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[Jh . Abstract 

We study the cosmological evolutions of the equation of state (EoS) for the universe in the ho- 
mogeneous and isotropic Priedmann-Leniaitre-Robertson- Walker (FLRW) space-time. In particu- 
lar, we reconstruct the cyclic universes by using the Weierstrass and Jacobian elliptic functions. 
R is explicitly illustrated that in several models the universe always stays in the non-phantom 






ff^ ' (quintessence) phase, whereas there also exist models in which the crossing of the phantom divide 

>■ 

can be realized in the reconstructed cyclic universes. 
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I. INTRODUCTION 



Recent observations on cosmic microwave background (CMB) radiation [l| have suggested 
that in the early universe inflation occurred and the universe became homogeneous and 
isotropic. On the other hand, according to various recent cosmological observations such as 



Type la Supernovae 



2|, CMB radiation 



IL large scale structure (LSS) [3|], baryon acoustic 
oscillations (BAO) j^, and weak lensing 5|, the current cosmic expansion is accelerating. 
Hence, in the present time so-called dark energy, which may be some unknown matter 
or geometrical energy coming from the deviation of a gravitational theory from general 
relativity, is dominant over non-relativistic matter, i.e., cold dark matter, radiation, and 
baryon (for reviews on dark energy, see, e.g., |6l-ll2|: for those on modifications of gravity. 



see, e.g. 



13 



Hl9|) 



In the beginning of the universe before inflation, it is considered that there was so-called 
a Big Bang singularity. In other words, the universe might have evolved from a singular 
space-time point. On the other hand, it is known that at dark energy dominant stage, there 
eventually appear the finite-time future singularities 20l-|25| . or there are several scenarios 
in which the late time universe contracts and eventually a Big Crunch singularity occurs. 
Thus, the issue discussed here is that there exist singularities in the beginning of the universe 
as well as at the last stage of it. If the evolution of the universe is periodic, the existence 
of a Big Bang singularity as well as the finite-time future singularities or a Big Crunch 
can be avoided. Based on this idea and inspired by string theories, the cyclic universe has 
been proposed in Ref. 26|]^ (for recent comprehensive papers on the cyclic universe, see, 
e.g., 28l-l3l|). In addition, the ekpyrotic scenario in the framework of the brane world has 
also been suggested in Ref. 32|. Furthermore, the bouncing universe has been investigated in 



Refs. 
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34l | (for a review on bouncing cosmology, see 35] ). Moreover, in Ref. 36| the (trefoil 



and figure-eight) knot universe has been studied, where the knot theory relates to the cyclic 
universe. In other words, the geometrical picture of the knot relates to oscillatory solutions 
of the gravitational field equations in the homogeneous and isotropic Friedmann-Lemaitre- 
Robert son- Walker (FLRW) and Bianchi-type I universes. Recently, there has also executed 
an investigation of figure eight knot in Ref. J37|. Furthermore, it is suggested 
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39| that 



^ In other context, the cychc universe has been argued in Ref. [2 



the Weierstrass p{t), ({t) and o"(t)-functions and the Jacobian elhptic functions can play 
an important role to examine astrophysical and cosmological problems (for recent studies 
of the applications, see, e.g., 40[). In particular, in Ref. |39[ the elliptic functions have 
been applied to describe the FLRW universe. Related cosmological features such as a cyclic 
behavior of the cosmic evolution have been studied in various scenarios 
Ref. 4^, the behavior of the EoS for dark energy has been investigated in the so-called 



4l|. Moreover, in 



g-essence models constructed by both k-essence 43|] and f-essence 4^, which is treated as 
a spinor field and corresponds to a classical c-number quantity. More recently, f-essence 
is dealt with a Grassmann- valued quantity in Ref. 45|] or an operator, i.e., q- number in 



Ref. 



46l | . Generalization of the Chaplygin gas ty pe models 47l | with the periodicity or the 



quasi-periodicity have also been explored in Ref. 48 1 



In this paper, with the Weierstrass and Jacobian elliptic functions, we reconstruct the 
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50(1 . It is important 



cyclic universe by using the equivalent procedure in Refs. [13|, 
to emphasize that to use the Weierstrass and Jacobian elliptic functions for describing the 
scale factor or the Hubble parameter is one of the novel ingredients in this work. The 
cosmological motivation to use such elliptic functions for the scale factor or the Hubble 
parameter is to realize a cyclic behavior of the universe naturally without special setting 
by using the properties of periodicities of the Weierstrass and Jacobian elliptic functions. 
This corresponds to the reconstruction procedures, through which an arbitrary cosmological 
expansion history can be reconstructed by starting with providing an appropriate form of 
the scale factor or the Hubble parameter with desirable features. We use the units of the 
gravitational constant SirG = c = 1 with G and c being the gravitational constant and the 
seed of light. 

The paper is organized as follows. In Sec. II, we explain the basic equations in the FLRW 
space-time. In Sec. HI, we explore models induced by the Weierstrass C(t)-function. In Sec. 
IV, we examine models induced by the Weierstrass cr(t)-function. In Sec. V, we investigate 
models induced by Jacobian elliptic functions. Finally, conclusions are given in Sec. VI. 



II. FLRW COSMOLOGY 

We start with the standard gravitational action 

S= I ^f^d^x{R + L^) , (2.1) 



where R is the scalar curvature and L^ is the Lagrangian of matter. Now, we assume the 
FLRW space-time with the metric, 

dr"^ 



ds^ = -dr + a\t) 



1-Kr 



2 



r'd^' 



(2.2) 



where a{t) is the scale factor, dQ'^ is the metric of 2-dimensional sphere with unit radius. 
Moreover, K is the cosmic curvature, and for K = —1, K = 0, and K = +1^, the universe is 
open, flat, and closed, respectively. Moreover, the Hubble parameter is defined hj H = a/ a, 
where a dot denotes the time derivative, d/dt. In the flat FLRW background (12. 2p . that is, 
k = 0, from the action in Eq. (12.11) . the gravitational equations and the continuity equation 
can be written in the H-form 

p = 3H\ (2.3) 

P = -2H - •iH'^ , (2.4) 

p=-3i7(p + P), (2.5) 

where p and P are the energy density and pressure, respectively. 

III. MODELS INDUCED BY THE WEIERSTRASS C(i)-FUNCTION 

In this section, we study a model (MG-II) induced by the Weierstrass C(^)-function. In 
particular, the Weierstrass (^(t)-function is the logarithmic derivative of the cr(t)-function, 
which we use in Sec. IV. We use the following procedure for the reconstruction of the 
EoS. First, we express the scale factor or the Hubble parameter by using the Weierstrass 
C(t)-function. Next, by combining these expressions with the gravitational field equations 
(12. 4 p and (12. Sp . we derive the energy density and pressure, and thus reconstruct the EoS. 
In this paper, the so-called Myrzakulov Gas (MG)-i {i = I, II, . . . , XV) model means 



K can take any value, but it is related to (— , 0, +) curvatures according to its sign. 



the model of gases or fluid, according to the notations of Ref. [36|]. We note that the 
Weierstrass C(t)-function is a quasi-periodic function and the Weierstrass p(t)-function is a 
two periodic function. Thus, in comparison with the use of ordinary trigonometric functions, 
the advantage of the use of the Weierstrass p(t)-function is that two periodic behaviors can 
be obtained. 

The physical motivation why we examine the MG-i gas is that in these models the cos- 
mological evolution of the generalized Chaplygin gas type models with the periodical and 
quasi-periodical features can be realized. The detailed behaviors are dependent on the 
models. The important point is that these models described by the Weierstrass functions. 
Therefore, the cosmic expansion history showing the periodicity and/or quasi-periodicity be- 
haviors can be demonstrated. As a consequence, these models can produce new cosmological 
scenarios in which cosmological singularities such as a Big Bang singularity, the finite-time 
future singularities or a Big Crunch singularity can be removed, similarly to that in, e.g., 
the cyclic universe, ekpyrotic and bouncing universe scenarios. 

A. MG-II model 

We investigate the MG-II model when the scale factor is written with the (^(t)-function 
as follows 

a = C(t) , (3.1) 

where ({t) is the (^-Weierstrass function. Hence, the Hubble parameter takes the form 

H-§. (3.) 

where p{t) is the p- Weierstrass function. In this case, the parametric EoS is given by 

P = 3(|||r, (3.3) 

where a prime denotes the time derivative of d/dt. The EoS parameter is expressed by 




FIG. 1: The scale factor a in Eq. (j3.ip as a function of t in the MG-II model for the Weierstrass 
invariants oi gi = 1 and §2 = 1. 



In Fig. [H we show the cosmological evolution of the scale factor a{t) in Eq. (13. ip as a 
function of t. We also depict the cosmological evolutions of the energy density p in Eq. fl3.3p 
and pressure P in Eq. (13. 4p as functions oft in Fig.O Furthermore, in Fig.[3]we demonstrate 
the cosmological evolution of the EoS w in Eq. fl3.5p as a function of t. Here, we have used 
the Weierstrass invariants oi gi = 1 and (72 = 1, which are defined to satisfy the following 
equation [39 1 

p-\t-gi,g2)= I ^dt'. (3.6) 



A{t'y-git'-g2 

From Fig. [H we see the oscillatory behavior of a. Thus, it is interpreted that by using 
the Weierstrass C(^)-function, a model describes the cyclic universe with two periods. This 
comes from the feature of the two periodicity of the Weierstrass (^(t)-function. We note that 
in Figs. [H and Ej there are no diverging behavior of a, p and p, namely, all the curves in each 
figure are smoothly connected^. 



On the other hand, the effective EoS for the universe in the FLRW space-time (2.2) is de- 
scribed by [isl Weff = Pcff/Pcff = -1-2H/ (3if2) with peff = 3i/2 and P^s = - (2H + 3hA, 



^ Since the amplitude is very large, in Figs. [T] and [2] apparently it seems there arc some divergence. Indeed, 
however, there are no divergence. 
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FIG. 2: The energy density p in Eq. (j3.3p [left panel] and pressure P in Eq. ([3 
functions of t. Legend is the same as Fig. [TJ 



[right panel] as 




FIG. 3: The EoS w in Eq. (j3.5p as a function of t. Legend is the same as Fig. [TJ 

where peff and PcS correspond to the total energy density and pressure of the universe, re- 
spectively. Throughout this paper, pcff is p in Eq. (2.2) and Pcff is P in Eq. (2.3). For H < 0, 
Wcs > — Ij which is the non-phantom (quintessence) phase, whereas for H > 0, Wcs < — 1, 
which is the phantom phase. For H = 0, Wcs = — 1, which corresponds to the cosmological 
constant. It is significant to remark that recent various observational data SJJ implies that 
the crossing of the phantom divide line of wbe = —1 occurred in the near past. Here, wbe 
is the EoS for dark energy and at the dark energy dominated stage, it can be regarded that 
w ~ i/^DE ~ Wcs- From Fig. [3], we find that multiple crossings of the phantom divide can be 
realized. 

We note that if the two periods of the elliptic functions are equal to infinity, that is, in 
the case rrii = m2 = oo, the elliptic functions are reduced (degenerate) to the elementary 



rational functions: 

a{t)=t, at)=r\ p{t)=t-\ (3.7) 

In this degenerate case, the MG-II model is reduced to the case 

a = -H = - (3.8) 

At the same time, the parametric EoS takes the form 

P = ^, P=-|. (3.9) 



This means that the EoS reads 



P = -^. (3.10) 



3 
Hence, the corresponding EoS parameter is given by 

w = --^-1.67, (3.11) 

so that the crossing of the phantom divide with w ~ —1.67 can be realized. Thus, we can 
conclude that the MG-II model is the two-periodical analogue of the usual universe filled by 
the barotropic fluid with the EoS of the form fl3.10p and that the crossing of the phantom 
divide with (13. lip occurs. 

IV. MODELS INDUCED BY THE WEIERSTRASS cT(t)-FUNCTION 

In this section, we explore two models (MG-I, MG-III) induced by the Weierstrass cr{t)- 
function, which is a quasi-periodic function. We reconstruct the EoS by using the same 
procedure as that in Sec. III. 

A. MG-I model 

We represent the scale factor as 

a = aocrit) , (4.1) 

where Oq is a constant and a(t) is the a- Weierstrass function. In this case, the Hubble 
parameter is described by 

H{t) = at) . (4.2) 
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FIG. 4: The EoS w in Eq. (j4.5p as a function of t. Legend is the same as Fig. [TJ 



Then, the gravitational field equations lead to the parametric EoS as 



P = 2p{t) - X\t) 



and the EoS parameter is given by 



w 



2pit)_ 



(4.3) 
(4.4) 



(4.5) 



In Fig. m we display the cosmological evolution of the EoS w in Eq. (14. 5 p as a function 
of t. Here, we have used the Weierstrass invariants oi gi = 1 and g2 = 1 in Eq. (13. 6p . It 
follows from Fig. HJthat the universe always stays in the non-phantom phase {w > —1). 

We also examine the case that two periods of the elliptic functions are equal to infinity, 
that is, in the case mi = m2 = oo, the elliptic functions are reduced (degenerate) to the 
elementary rational functions according to the equations in (13. 7p . In this degenerate case, 
the MG-I model is reduced to the case in which the scale factor and the Hubble parameter 
are given by 

a = a^t , H = - . (4.6) 



At the same time, the parametric EoS is expressed as 

3 „ 1 



P = ^, P 



t^' 



(4.7) 



Therefore, the EoS is described as 



^ = -f- (4.8) 



Hence, the corresponding EoS parameter becomes 

^ = '1- (4-9) 

As a result, it can be concluded that the MG-I model is the two-periodic analogue of the 
usual universe filled with the barotropic fluid with the EoS in Eq. f l4.8p and with the EoS 
parameter in Eq. fl4.9p . 

B. MG-III model 

Next, we express the Hubble parameter as 

H{t) = (j{t) . (4.10) 

From this representation of H, the scale factor is given by 

a = aoexp( I a{K)dKj . (4.11) 

Thus, the gravitational field equations give the parametric EoS as 

p = 3a\t) , (4.12) 

P = ~2a{t)C{t) - 3a\t) , (4.13) 

Moreover, the EoS parameter becomes 

» = -l-HM. (4.14) 

3 a(t) ^ ' 

In Fig. [5l we show the cosmological evolutions of the scale factor a in Eq. (14. lip and the EoS 
w in Eq. (I4.14P as functions of t for the Weierstrass invariants oi gi = 1 and g2 = I- From 
Fig. El we see that there does not exist the oscillating behavior of a, and that crossings of 
the phantom divide can be realized. 

Now, we consider the case that two periods of the elliptic functions are equal to infinity, 
i.e., nil = '^2 = oo, the elhptic functions are reduced (degenerate) to the elementary rational 
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FIG. 5: The scale factor a in Eq. (j4.1ip [left panel] and the EoS w in Eq. (j4.14p [right panel] as 
functions of t in the MG-III model. Legend is the same as Fig. [TJ 



functions in (13. 7p . In this degenerate case, the MG-III model is reduced to the model in 
which the Hubble parameter and the scale factor are given by 

H = t, a = aoe°-^*' . 



Moreover, the parametric EoS becomes 

p = 3t2 , P = -2-3f. 
It follows from these equations that the EoS is written as 

P = -2- p. 
Therefore, the corresponding EoS parameter reads 



w 



'^ 3t2 ■ 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



Consequently, Eq. ( J4.18P informs us that the EoS parameter is always less than —1 as 
w < —1 and in the limit of t ^ oo, w = —1. Thus, the late-time accelerated expansion of 
the universe can be realized. 



V. MODELS INDUCED BY THE JACOBIAN ELLIPTIC FUNCTIONS 

In this section, we investigate FLRW models (MG-V, MG-VI, . . . , MG-XVI) induced by 
the Jacobian elliptic functions cnt = cn(t, ?7i), snt = sn(t, m) and dnt = dn(t, ?7i), where m 
is the parameter of the elliptic modulus. The Jacobian elliptic functions are doubly periodic 
functions. Hence, these functions lead to a new class of cosmological models of the cyclic 
universes with a two periodic feature. With the same procedure as that in Sees. HI and IV, 
we reconstruct the EoS. 
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A. Formulations 



First, we present some differential equations for the Jacobian elliptic functions 52| 



sn't = cntdnt , 

cn't = — sntdnt , 

dn't = —rn^sntcnt , 

sn"t = — sntdn^t — TTi^sntcn^t = — snt(dn^t + m^cn^t) 

cn"t = cnt(?Ti^sn^t — dn t) , 

dn"t = m^dnt(sn^t — cn^t) , 



(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 



where cn't = dcnt/dt and so on. Hence, we get some known equations for these functions. 
For example, the differential equations for y = snt reads 52 1 



dt 



2 + (1 + m^)y - 2m^y^ 







^) =(l-l/^)(l-mV). 



For y = cnt, we have the differential equations 52 1 

d'y 



dt2 

dy 



1 - 2m^)y + 2m\^ = , 



^) =il-y')il-m' + mY). 



(5.7) 
(5.8) 

(5.9) 
(5.10) 



(2 - m^)y + 2y' 







In addition, we present the differential equations for the function y = dnt 52 1 

d^ 

dt2 

dy^ 

dt, 

We note that these functions are doubly periodic generalizations of the trigonometric func- 
tions and give them for the particular values of the parameter m: 



^] =(y^-l)(l-m'-y') 



(5.11) 
(5.12) 



sn(t, 0) = sin(t) , cn(t, 0) = cos(t) , dn(t, 0) = 1 . 
Also, we present other degenerate values of these functions: 

sn(t, 1) = tan(t) , cn(t, 1) = sec(t) , dn(t, 1) = sec(t) . 



(5.13) 



(5.14) 



In what follows, we demonstrate several examples of the cyclic cosmological models induced 
by the Jacobian elliptic functions. 
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FIG. 6: The scale factor a in Eq. (j5.15p as a function of t in the MG-V model for the elliptic 
modulus of m = 2. 



B. MG-V model 



As an example, for MG-V model, the scale factor is represented as 



a{t) = a^e 



cn(i) 



(5.15) 



The Hubble parameter is described by 



H = cn'(t) = -sn(t)dn(t) 



(5.16) 



With the gravitational field equations, the parametric EoS is expressed as 

p = 3(sn(t)dn(t))2 , 
P = -3dn^{t)sn^{t) + 2cn(t) (dn^it) - nsn^{t)) . 



(5.17) 
(5.18) 



The EoS parameter is given by 



w = -1 + -cn(t) 
o 

where we have used n-elliptic modulus. 



n 1 

+ 



dn\t) sn2(t) 



2{+\ y ' 



(5.19) 
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FIG. 7: The energy density p in Eq. (j5.17p [left panel] and pressure P in Eq. ()5.18|) [right panel] 
as functions of t. Legend is the same as Fig. [6l 
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FIG. 8: The EoS w in Eq. ()5.19p as a function of t. Legend is the same as Fig. [6l 



We demonstrate the cosmological evolution of the scale factor a{t) in Eq. fl5.15p as a 
function of t in Fig. [6l We also show the cosmological evolutions of the energy density p in 
Eq. fl5.17p and pressure P in Eq. flS.lSp as functions of t in Fig. [71 Furthermore, in Fig. [8] 
we illustrate the cosmological evolution of the EoS w in Eq. f l5.19p as a function of t. Here, 
we have taken the parameter of the elliptic modulus parameter as m = 2. From Fig. [HI we 
see the oscillation of a. Thus, it is considered that by using the Jacobian elliptic functions, 
a model with realizing the cyclic universe can be reconstructed. Moreover, it is clearly seen 
from Fig. [8] that crossings of the phantom divide can be realized. 

In Appendix A, for other FLRW models (MG-VI, . . . , MG-XVI), we present the scale 
factor, the energy density and pressure, and the parametric EoS as well as the EoS parameter. 
We have executed the same analysis as that shown in this section. In each model, we have 
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examined the cosinological evolutions of the scale factor and the EoS as functions of t. 

VI. CONCLUSIONS 

In the present paper, we have explored the cosmological evolutions of the EoS for the uni- 
verse in the homogeneous and isotropic FLRW background. With the Weierstrass p{t), ((t) 
and cr(t)-functions and the Jacobian elliptic functions, we have reconstructed cosmological 
models which can describe the cyclic universes. Furthermore, we have explicitly demon- 
strated that in the MG-I, MG-XI, MG-XXXIV, MG-XXXVI, MG-XXXVII, MG-XXXVIII 
models, the universe always stays in the non-phantom (quintessence) phase, whereas in the 
MG-II, MG-III, MG-V, MG-VI, MG-VII, MG-VIII, MG-IX, MG-X, MG-XXXV models, 
the crossing of the phantom divide can be realized. 

It is known that there exist two approaches to produce the cyclic universes. One is to 
introduce non-canonical scalar field, which makes the vacuum unstable in the quantum field 
theory. Another is to extend the gravitational theory, such as gravity with the higher order 
derivative term and f{R) gravity. It is also important to indicate that our MG-i (i=I, II, 
■ ■ ■ ) models are 2-periodic or 1-periodic generalizations of the usual non-periodic FLRW 
models. 

We note that the issue of removing singularities is a fundamental question. In other words, 
this is what a new physics comes out, when in the cosmological models the energy/ curvature 
scale is so high that render the hypothesis of general relativity should be rendered inapplica- 
ble. Although this is a question far from a satisfactory solution, a complete picture of these 
cosmologies certainly will need to take this issue seriously into account. Furthermore, there 
exists models, in particular, based on higher-curvature theories, with possessing bouncing 
solutions, i.e., without curvature singularities, that modify general relativity only if the 
universe becomes close to the predicted singularity, e.g., at the Planck scale. On the one 
hand, these theories would pass the solar test equally well as general relativity. Moreover, 
in principle, the perspective of cyclic cosmologies could be realized without supposing the 
existence of singularities as an unavoidable conclusion of the theory. This issue should be 
considered in more detailed in the future works. 



In Ref. 12|, it has been verified that all the dark energy cosmologies can be realized by 



different fluids of the universe with a general form of the EoS, and that the evolutions of all 
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the dark energy universes at the present time can be similar to that of the ACDM model, 
in which the universe consists of the cosmological constant A and cold dark matterfCDM). 
Since the ACDM model is compatible with the various cosmological observations iHSj, the 
models of the dark energy fluids are also consistent with the current observations. Indeed, it 
has been illustrated that different dark energy models are equivalent by examining various 
scalar field theories such as single and multiple scalar fields theories as well as tachyon scalar 
theory and holographic dark energy. In these theories, the current accelerated expansion 
of the universe with the quintessence or phantom phases can be realized. Also, the equiva- 
lence of these theories to the corresponding fluid descriptions have been shown. According 
to these consequences, we can understand that for the MG-I, MG-XI, MG-XXXIV, MG- 
XXXVI, MG-XXXVII, MG-XXXVIII models, the non-phantom (quintessence) phase occurs 
and hence these fluid models correspond to the quintessence models. On the other hand, for 
the MG-II, MG-III, MG-V, MG-VI, MG-VII, MG-VIII, MG-IX, MG-X, MG-XXXV models, 
since the crossing of the phantom divide can happen, these fluid models can be re gard ed to 
be equivalent to the two scalar fields models, for example, the quintom models J8|, |57| with 
both the canonical and non-canonical scalar fields. 

The new cosmological consequence obtained in this work is that by combining the re- 
construction method of the EoS for the universe with mathematical special functions such 
as the Weierstrass and Jacobian elliptic functions, new descriptions of the cyclic universes 
can be acquired. It is considered that this can be worthy of a clue of one of the generalized 
describing manner of the expansion history of the universe. 

In addition, it is remarkable to mention that if we apply the investigations in this work 
to the spatially anisotropic space-time, the EoS for the universe has the asymmetry during 
the oscillatory process of the expansion and contraction, so that the cosmological hysteresis 
can be realized 31 1. 



Finally, it is also important to stress that concerning the issue on potential deviation of 
a gravitational theory from general relativity, the detection of gravitational waves could be 
a definitive tool for discrimination between a gravitational theory and general relativity, as 



shown in Ref. 58 1. 
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Appendix A: Other models 

In this appendix, we describe other FLRW models (MG-VI, ..., MG-XVI). For each 
model, by executing the numerical calculations for the parameter of the elliptic modulus of 
r/i = 2, we have analyzed the cosmological evolutions of the scale factor and the EoS. 

1. MG-VI model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 

a{t) = aoe^^W , (Al) 

H = sn'(t) = dn(t)cn(t) , (A2) 

p = 3[cn(t)dn(t)]^ (A3) 

P = -3dn\t)cn^{t) - 2sn(t) {-dn^{t) - ncn'^{t)) , (A4) 

We have confirmed the oscillation of a and found that crossings of the phantom divide can 
be reahzed. 
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2. MG-VII model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 

a(t) = aoe^^W , (A6) 

H = dn't , (A7) 

p = 3[ncn(t)sn(t)]2 , (A8) 

P = -3n\n^{t)cn^{t) + 2ndn(t) {cn^{t) - sn^{t)) , (A9) 

2d^/ 1 _ 1 y 

We have confirmed the oscillation of a and found that crossings of the phantom divide can 
be realized. 

3. MG-VIII model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 

a{t) = aoe.p(^^^^^^^^] , (All) 

V ^i-dn^(t) ; 

H = cnt, (A12) 

p = 3cn2(t), (A13) 

P = -3cn2(t) + 2dn(t)sn(t) , (A14) 

.._l,?^^(^. ,A15) 

3cn^(t) 

We have confirmed that there is no oscillating behavior of a and found that crossings of the 
phantom divide can be realized. 
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4. MG-IX model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 

a{t) = ao(-v^cn(t) + dn(t))^/v^ . (A16) 

H = sn(t) , (A17) 

p = 3sn2(t) , (A18) 

P = -3sn2(t) - 2dn(t)cn(t) , (A19) 

We have seen that crossings of the phantom divide can be reahzed. 

5. MG-X model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 

a{t) = aoe"'"W , (A21) 

H = dnt, (A22) 

p = 3dn^{t) , (A23) 

P = -3dn2(t) + 2nsn(t)cn(t) , (A24) 

We have confirmed the oscillation of a and found that crossings of the phantom divide can 
be realized. 
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6. MG-XXXIII model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 



a{t) = cn'(t) = -dn(t)sn(t) . 
/dn(t) sn(t)\ 



H = cnt 



Vsn(t) dn(t)y 



p^3(cntf^-n^y^ 
' \sn{t) dnit)J, 

^'\sn{t) dn{t)J 



(A26) 
(A27) 

(A28) 



P 



-dn\t)sn\t) + ndn\t)sn\t) - cn\t) {dn\t) + nsn'^it))' 

dn2(t)sn2(t) 

-dn'*(t)sn2 (t)+ndn^ (t)sn4 (t)-cn2 (t) ("dn^ (t)+nsn2 (t)) 
dn^(t)sn2(t) 



w 



dn{t) ^ sn(f) 



(A29) 



(A30) 



We have confirmed the oscillation of a and found that the universe always stays in the 
non-phantom (quintessence) phase {w > —1). 
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7. MG-XXXIV model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 



a{t) = sn'(t) = cn(t)dn(t) , 

sn(t) (racn2(t) + dn^(t)) 

cn(t)dn(t) 

^sn2(t) (ncn2(t) + dn^(t))' 



H 



P 



cn2(t)dn'(t) 
,sn2(t)(ncn2(t) + dn2(t))' 
' cn2(t)dn2(t) 

-dn^{t) + sn^{t) (271- 



dn\t) 
cn2(t) 



w 



-dn\t) + sn\t) l^2n - ^^ 



ncn^(t) j — 1 — n 
+ ncn^{t) \ —1 — n 



sn2(t) 



dn\t) 



dn (t) 



^ ^ dn (t)cn2(t) 



(A31) 
(A32) 

(A33) 



(A34) 



(A35) 



We have confirmed the oscillation of a and found that the universe always stays in the 
non-phantom (quintessence) phase {w > —1). 

8. MG-XXXV model 



The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 



a{t) = —ncn(t)sn(t) , 

cn(t)dn(t) sn(t)dn(t) 



H 



p = 3 



sn(t) 



cn(t) 



cn(t)dn(t) sn(t)dn(t) 



sn(t) 



cn(t) 



cn(t)dn(t) sn(t)dn(t) 



sn(t) 



cn(t) 



-2dn^(t) + cn^(t) ( -n - ^) + sn\t) (n ^"'^'^ 



w 



2dn^(t) + cn\t) ^-n - ^) + sn\t) (n - ^ 



o / cn(t)dn(t) _ sn(f)dn(t) 
■^ ' sn(t) cn(t) 



(A36) 
(A37) 

(A38) 



(A39) 



(A40) 
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We have confirmed the oscillation of a and found that the crossings of the phantom divide 
can be realized. 

9. MG-XXXVI model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 



a{t) = cn(t) , 

sn(t)dn(t) 



H 



p = 3 



P 



w 



cn(t) ' 

sn(t)dn(t) 

cn(t) 

sn(t)dn(t) 
cn(t) 



nsn^it) + dn^(t) 



sn2(t) 
cn2(t) 



nsn^it) + dn^(t) 



sn^(t) 
cn2(t) 



cn2(t) 



q / sn(t)dn(f) 
"^ ' cn(t) 



(A41) 
(A42) 

(A43) 

(A44) 

(A45) 



We have confirmed the oscillation of a and found that the universe always stays in the 
non-phantom (quintessence) phase {w > —1). 

10. MG-XXXVII model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 



a(t) = sn(t) , 

cn(t)dn(t) 



H 



p = 3 



P = -3 



sn(t) ' 
cn(t)dn(t) 
sn(t) 

cn(t)dn(t) 
sn(t) 



-dn\t) + cn^{t) 



-n 



dn\t) 
sn2(t) 



w 



-1 - 



-dn^(t) + cn2(t) 



-n — 



dn^(f) 

sn2(t) 



^ cn(f)dn(t) ' 
sn(i) 



(A46) 
(A47) 

(A48) 

(A49) 

(A50) 
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We have confirmed the oscillation of a and found that the universe always stays in the 
non-phantom (quintessence) phase {w > —1). 

11. MG-XXXVIII model 

The scale factor, the energy density and pressure, and the parametric EoS, and the EoS 
parameter are given by 



a(t) = dn(t) , 

cn(t)sn(t) 



H = -n- 



p = 3 [n 



dn(t) ' 

cn(t)sn(t) 

dn(t) 



P 



w 



_ cn(t)sn(t)| _2^ 



dn(t) 



sn2(t) +cn^{t) i-l-n 



sn2(t) 
'd^ 



2n 



sn2(t) + cn2(t)(-l-n|^) 



o /^ cn(f)sn(t) y 
■^"^ dnw ) 



(A51) 
(A52) 

(A53) 

(A54) 

(A55) 



We have confirmed the oscillation of a and found that the universe always stays in the 
non-phantom (quintessence) phase {w > —1). 



[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) 



|arXiv:astro-ph/0302209 1 ; Astrophys. J. Suppl. 170, 377 (2007) |arXiv:astro-ph/0603449| ; 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) 

|arXiv:0803.0"547] [astro-ph]]; Astrophys. J. Suppl. 192, 18 (2011) |arXiv:1001.4538l [astro- 

ph.CO]]; G. Hinshaw, D. Larson, E. Komatsu, D. N. Spergel, C. L. Bennett, J. Dunkley, 

M. R. Nolta and M. Halpern et a/. , |ar2<iYil212,5226j [astro-ph. CO]. 
[2] S. Perlmutter et al. [SNCP Collaboration], Astrophys. J. 517, 565 (1999) 

[arXiv:astro-ph/9812133] ; A. G. Riess et al. [Supernova Search Team Collaboration], 

Astron. J. 116, 1009 (1998) |arXiv:astro-ph/9805201] . 
[3] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) 

|arXiv:astro-ph/0310723] ; U. Seljak et al. [SDSS Collaboration], Phys. Rev. D 71, 103515 

(2005) |arXiv:astro-ph/0407372] . 



23 



[4] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633, 560 (2005) 
|arXiv:astro-ph/050lT7l] . 

[5] 



[6] 

[7] 

[s; 

[9 
[10 

[11 

[12 

[13 

[14 
[15 
[16 
[17 
[18 

[19 
[20 



B. Jain and A. Taylor, Phys. Rev. Lett. 91, 141302 (2003) [arXiv:astro-ph/0306046] . 

E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006) 

|arXiv:hep-th/0603057] . 

R. Durrer and R. Maartens, Gen. Rel. Grav. 40, 301 (2008) |arXiv:071 1.0077) [astro-ph]]; 

larXiv:0811.4T32l l [astro-ph]. 

Y. F. Cai, E. N. Saridakis, M. R. Setare and J. Q. Xia, Phys. Rept. 493, 1 (2010) 

|arXiv:0909.27761 [hep-th]] . 

S. Tsuiikawa. larXiv: 1004. 1493] [astro-ph.CO]. 

L. Amendola and S. Tsujikawa, Dark Energy (Cambridge University press, 2010). 

M. Li, X. D. Li, S. Wang and Y. Wang, Commun. Theor. Phys. 56, 525 (2011) [arXiv:1103.5870l 

[astro-ph.CO]]. 

K. Bamba, S. Capozziello, S. Nojiri and S. D. Odintsov, Astrophys. Space Sci. 342, 155 (2012) 



|arXiv:1205.342ll [gr-qc]]. 

S. Nojiri and S. D. Odintsov, Phys. Rept. 505, 59 (2011) [arXiv:1011.0"544] [gr-qc]]; eConf 

C0602061, 06 (2006) [Lit. J. Geom. Meth. Mod. Phys. 4, 115 (2007)] |arXiv:hep-th/0601213] . 

T. P. Sotiriou and V. Faraoni, Rev. Mod. Phys. 82, 451 (2010) [arXiv:0805. 17261 [gr-qc]]. 

S. Capozziello and V. Faraoni, Beyond Einstein Gravity (Springer, 2010). 

S. Capozziello and M. De Laurentis, Phys. Rept. 509, 167 (2011) [arXiv: 1108.6266) [gr-ac]]. 

A. De Felice and S. Tsujikawa, Living Rev. Rel. 13, 3 (2010) [arXiv: 1002.4928) [gr-ac]]. 

T. Clifton, P. G. Ferreira, A. Padilla and C. Skordis, Phys. Rept. 513, 1 (2012) 

|arXiv:1106.2476) [astro-ph.CO]]. 

S. Capozziello, M. De Laurentis and S. D. Odintsov. larXiv: 1206.484? [gr-qc]. 

R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 071301 

(2003) [arXiv:astro-ph/0302506| ; B. Mclnnes, JHEP 0208 (2002) 029 |arXiv:hep-th/0112066] ; 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 562, 147 (2003) |arXiv:hep-th/0303117] ; 

Phys. Lett. B 571, 1 (2003) |arXiv:hep-th/0306212j ; V. Faraoni, Int. J. Mod. Phys. D 11, 



471 (2002) |arXiv:astro-ph/0110067 J; P. F. Gonzalez-Diaz, Phys. Lett. B 586, 1 (2004) 
|arXiv:astro-ph/0312579l ; E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 
043539 (2004) |arXiv:hep-th/0405034| ; P. Singh, M. Sami and N. Dadhich, Phys. Rev. D 



24 



68, 023522 (2003) [arXiv:hep-th/0305110| ; C. Csaki, N. Kaloper and J. Terning, Annals 
Phys. 317, 410 (2005) |arXiv:astro-ph/0409596| ; P. X. Wu and H. W. Yu, Nucl. Phys. 
B 727, 355 (2005) [arXiv:astro-ph/0407424] ; S. Nesseris and L. Perivolaropoulos, Phys. 
Rev. D 70, 123529 (2004) |arXiv:astro-ph/0410309l ; M. Sami and A. Toporensky, Mod. 
Phys. Lett. A 19, 1509 (2004) [arXiv:gr-qc/0312009] ; H. Stefancic, Phys. Lett. B 586, 
5 (2004) |arXiv:astro-ph/0310904] ; L. P. Chimento and R. Lazkoz, Phys. Rev. Lett. 91, 
211301 (2003) |arXiv:gr-qc/0307111J ; J. G. Hao and X. Z. Li, Phys. Lett. B 606, 7 (2005) 
|arXiv:astro-ph/0404154] ; E. Ehzalde, S. Nojhi, S. D. Odintsov and P. Wang, Phys. Rev. 



D 71, 103504 (2005) |arXiv:hep-th/0502082| ; M. P. Dabrowski and T. Stachowiak, Annals 
Phys. 321, 771 (2006) |arXiv:hep-th/0411199] ; F. S. N. Lobo, Phys. Rev. D 71, 084011 



(2005) |arXiv:gr-qc/0502099i ; R. G. Cai, H. S. Zhang and A. Wang, Commun. Theor. Phys. 
44, 948 (2005) |arXiv:hep-th/0505186] ; I. Y. Aref'eva, A. S. Koshelev and S. Y. Vernov, 
Phys. Rev. D 72, 064017 (2005) |arXiv:astro-ph/0507067| ; H. Q. Lu, Z. G. Huang and 



W. Fang, |arXiv:hep-th/0504038[ W. Godlowski and M. Szydlowski, Phys. Lett. B 623, 10 
(2005) [arXiv:astro-ph/0507322| ; J. Sola and H. Stefancic, Phys. Lett. B 624, 147 (2005) 
|arXiv:astro-ph/0505133] ; B. Guberina, R. Horvat and H. Nikohc, Phys. Rev. D 72, 125011 
(2005) |arXiv:astro-ph/0507666| . 



[21] Y. Shtanov and V. Sahni, Class. Quant. Grav. 19, LlOl (2002) |arXiv:gr-qc/0204040| . 

[22] J. D. Barrow, Class. Quant. Grav. 21, L79 (2004) |arXiv:gr-qc/0403084| . 

[23] S. Nojiri and S. D. Odintsov, Phys. Lett. B 595, 1 (2004) |arXiv:hep-th/0405078] ; 
Phys. Rev. D 70, 103522 (2004) |arXiv:hep-th/0408170| ; Phys. Rev. D 72, 023003 
(2005) [arXiv:hep-th/0505215| ; S. Cotsakis and I. Klaoudatou, J. Geom. Phys. 55, 
306 (2005) |arXiv:gr-qc/040902 2|; M. P. Dabrowski, Phys. Rev. D 71, 103505 (2005) 
|arXiv:gr-qc/0410033 |; L. Fernandez- Jambrina and R. Lazkoz, Phys. Rev. D 70, 121503 (2004) 
|arXiv:gr-qc/0410124| ; Phys. Lett. B 670, 254 (2009) |arXiv:0805. 22841 [gr-ac]]: J. D. Barrow 
and C. G. Tsagas, Class. Quant. Grav. 22, 1563 (2005) |arXiv:gr-qc/0411045] ; H. Stefancic, 
Phys. Rev. D 71, 084024 (2005) |arXiv:astro-ph/0411630 | ; C. Cattoen and M. Visser, Class. 
Quant. Grav. 22, 4913 (2005) |arXiv:gr-qc/0508045| ; P. Tretyakov, A. Toporensky, Y. Shtanov 



and V. Sahni, Class. Quant. Grav. 23, 3259 (2006) (arXiv:gr-qc/0510104j ; A. Balcerzak and 



M. P. Dabrowski, Phys. Rev. D 73, 101301 (2006) |arXiv:hep-th/0604034] ; M. Sami, P. Singh 
and S. Tsujikawa, Phys. Rev. D 74, 043514 (2006) |arXiv:gr-qc/0605113| ; M. Bouhmadi-Lopez, 



25 



p. F. Gonzalez-Diaz and P. Martin-Moruno, Phys. Lett. B 659, 1 (2008) |arXiv:gr-qc/0612135] ; 

A. V. Yurov, A. V. Astashenok and P. F. Gonzalez-Diaz, Grav. Cosmol. 14, 205 (2008) 

|arXiv:0705.4T08] [astro-ph]] : T. Koivisto, Phys. Rev. D 77, 123513 (2008) [arXiv:0803.3"399] [gr- 

qc]]; J. D. Barrow and S. Z. W. Lip, Phys. Rev. D 80, 043518 (2009) [ arXiv:0901. 16261 [gr-ac]]: 

M. Bouhmadi-Lopez, Y. Tavakoli and P. V. Moniz, JCAP 1004, 016 (2010) ; arXiv:0911.i428] 

[gr-qc]]. 
[24] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71 (2005) 063004 

|arXiv:hep-th/05010"25] . 
[25] S. Nojiri and S. D. Odintsov, Phys. Rev. D 78, 046006 (2008) |arXiv:0804.3"5T9] [hep-th]]; 

K. Bamba, S. Nojiri and S. D. Odintsov, JCAP 0810 (2008) 045 [arXiv:0807.2"575] [hep-th]]; 

K. Bamba, S. D. Odintsov, L. Sebastiani and S. Zerbini, Eur. Phys. J. C 67 (2010) 295 

|arXiv:0911.4"390] [hep-th]]. 
[26] P. J. Steinhardt and N. Turok, Phys. Rev. D 65, 126003 (2002) |arXiv:hep-th/01 11098] ; 



arXiv:hep-th/01 11030 J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. 



D 66, 046005 (2002) |arXiv:hep-th/ 0109050|; P. J. Steinhardt and N. Turok, Nucl. Phys. 
Proc. Suppl. 124, 38 (2003) |arXiv:astro-ph/0204479| ; Science 296, 1436 (2002); J. Khoury, 
P. J. Steinhardt and N. Turok, Phys. Rev. Lett. 92, 031302 (2004) |arXiv:hep-th/0307132] ; 
P. J. Steinhardt and N. Turok, Science 312, 1180 (2006) |astro-ph/0605173] . 



[27] D. Y. Chung, arXiv:physics/0105064 



[28] K. Saaidi, H. Sheikhahmadi and A. H. Mohammadi, Astrophys. Space Sci. 338, 355 (2012) 
|arXiv:1201. 0275] [gr-ac]]: S. Nojiri, S. D. Odintsov and D. Saez-Gomez, AIP Conf. Proc. 1458, 
207 (2011) |arXiv:1108.0767l [hep-th]]. 

[29] Y. F. Cai and E. N. Saridakis, J. Cosmol. 17, 7238 (2011) [ arXiv:1108.6052l [gr-ac]]. 

[30] T. Biswas and A. Mazumdar, Phys. Rev. D 80, 023519 (2009) |arXiv:0901.4930] [hep-th]]; 
T. Biswas, T. Koivisto and A. Mazumdar, JCAP 1011, 008 (2010) [arXiv:1005.0590] [hep-th]]: 
T. Biswas, E. Gerwick, T. Koivisto and A. Mazumdar, Phys. Rev. Lett. 108, 031101 (2012) 
|arXiv:1110.5"249] [gr-qc]]. 

[31] V. Sahni and A. Toporensky, Phys. Rev. D 85, 123542 (2012) [arXiv: 1203.03951 [gr-ac]]. 

[32] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 (2001) 
|arXiv:hep-th/0103239] ; |arXiv:hep-th/0105212t R. Y. Donagi, J. Khoury, B. A. Ovrut, 



P. J. Steinhardt and N. Turok, JHEP 0111, 041 (2001) [arXiv:hep-th/0105199| ; J. Khoury, 

26 



B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 66, 046005 (2002) 
|arXiv:hep-th/01090"50] . 
[33] D. N. Page, Class. Quant. Grav. 1, 417 (1984). P. Peter and N. Pinto-Neto, Phys. Rev. D 65, 



023513 (2001) | arXiv:gr-qc/0109038| ; 66, 063509 (2002) |arXiv:hep-th/0203013| ; Y. Shtanov 
and V. Sahni, Phys. Lett. B 557, 1 (2003) |arXiv:gr-qc/0208047 ]; T. Biswas, A. Mazumdar 
and W. Siegel, JCAP 0603, 009 (2006) |arXiv:hep-th/0508194| ; Y. F. Cai, T. Qiu, Y. S. Piao, 
M. Li and X. Zhang, JHEP 0710, 071 (2007) |arXiv:0704. 10901 [gr-qc]]; P. Creminehi and 



L. Senatore, JCAP 0711, 010 (2007) |arXiv:hep-th/0702165| ; Y. S. Piao, Phys. Lett. B 677, 
1 (2009) [arXiv:0901.2"644 i [gr-qc]]: J. Zhang, Z. G. Liu and Y. S. Piao, Phys. Rev. D 82, 123505 
(2010) [arXiv: 1007.24981 [hep-th]]: Y. S. Piao, Phys. Lett. B 691, 225 (2010) [ arXiv:1001.063T] 
[hep-th]]; Z. -G. Liu and Y. -S. Piao, Phys. Rev. D 86, 083510 (2012) [arXiv:1201.i37T] [gr- 
qc]]; P. P. Avehno and R. Z. Ferreira, Phys. Rev. D 86, 041501 (2012) [arXiv: 1205.6676) 
[astro-ph.CO]]. 

[34] C. Corda and H. J. Mosquera Cuesta, Astropart. Phys. 34, 587 (2011) [ arXiv:1011.480T] 
[physics. gen-ph]] . 

[35] M. Novello and S. E. P. Berghaffa, Phys. Rept. 463, 127 (2008) jarXiv: 0802. 1634] [astro-ph]]. 

[36] R. Myrzakulov, Gen. Rel. Grav. 44, 3059 (2012) [arXiv: 1008.4486] [physics. gen-ph]]; 
larXiv:1204.T093i [phvsics. gen-ph]: R. Myrzakulov et al.. Int. J. Theor. Phys. 51, 1204 (2012) 
[arXiv: 1102.4456 1 [phvsics.gen-ph]]: [arXiv: 1201.43601 [physics. gen-ph]. 

[37] H. Itoyama, A. Mironov, A. Morozov and A. .Morozov, JHEP 1207, 131 (2012) 
|arXiv:1203.5978] [hep-th]]. 

[38] G. W. Gibbons and M. Vyska, Class. Quant. Grav. 29, 065016 (2012) [ arXiv: 11 10.6508] 
[gr-qc]]; L Bochicchio, S. Capozziello and E. Laserra, Lit. J. Geoni. Meth. Mod. Phys. 
8, 1653 (2011) far Xiv: 1111.0819 ^ [gr-qc]]; B. G. Dimitrov, J. Math. Phys. 44, 2542 (2003) 
|arXiv:hep-th/01072"3T] . 

[39] J. D'Ambroise,|ar2£iYiQ908,2481J [gr-qc]. 

[40] R. Teh, P. -Y. Tan and K. -M. Wong, Int. J. Mod. Phys. A 27, 1250148 (2012) jarXiv: 1203.03821 
[hep-th]]; V. P. Spiridonov and G. S. Vartanov, JHEP 1206, 016 (2012) [arXiv:1203.5677l [hep- 
th]]. 

[41] M. Jamil, N. A. Myrzakulov, K. K. Yerzhanov, D. Momeni and R. Myrzakulov, 
I arXiv: 1201.43601 [physics. gen-ph]; A. J. Lopez- Revelles, R. Myrzakulov and D. Saez-Gomez, 

27 



Phys. Rev. D 85, 103521 (2012) |arXiv:1201.5647] [gr-ac]]: M. Sharif, K. R. Yesmakhanova, 

S. Rani and R. Mvrzakulov. larXiv: 1204.2181 1 [phvsics.gen-ph]: R. Mvrzakulov. larXiv: 1205.52661 

[physics.gen-ph]. 
[42] K. Bamba, O. Razina, K. Yerzhanov and R. Mvrzakulov. larXiv:1203. 28041 [phvsics.gen-ph]. 
[43] T. Chiba, T. Okabe and M. Yamaguchi, Phys. Rev. D 62, 023511 (2000) 

[arXiv:astro-ph/9912463| ; C. Armendariz-Picon, T. Damour and V. F. Mukhanov, Phys. Lett. 

B 458, 209 (1999) |arXiv:hep-th/9904075] ; J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 



219 (1999) [ arXiv:hep-th/9904176] ; C. Armendariz-Picon, V. F. Mukhanov and P. J. Stein- 
hardt, Phys. Rev. Lett. 85, 4438 (2000) |arXiv:astro-ph/0004134] ; Phys. Rev. D 63, 103510 
(2001) |arXiv:astro-ph/0006373] ; R. de Putter and E. V. Linder, Astropart. Phys. 28, 263 
(2007) |arXiv:0705.0400l [astro-ph]]. 

[44] R. Mvrzakulov. larXiv:1011.4337] [astro-ph.CO]. 

[45] T. Damour and C. Hillmann, JHEP 0908, 100 (2009) [ arXiv:0906.3TT6] [heo-th]]. 

[46] T. Damour and P. Spindel, Phys. Rev. D 83, 123520 (2011) |arXiv:1103.2927 ' [gr-ac]]. 

[47] A. Y. Kamenshchik, U. Moscheha and V. Pasquier, Phys. Lett. B 511 (2001) 265 
|arXiv:gr-qc/0103004| ; M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev. D 66 (2002) 
043507 [arXiv:gr-qc/0202064]. 



[48] K. Bamba, U. Debnath, K. Yesmakhanova, P. Tsyba, G. Nugmanova and R. Myrzakulov, 

Entropy 14, 2351 (2012) |arXiv:1203.4"226] [gr-ac]]. 
[49] S. Nojiri and S. D. Odintsov, Phys. Rev. D 72 (2005) 023003 |arXiv:hep-th/0505215| . 
[50] H. Stefancic, Phys. Rev. D 71 (2005) 084024 |arXiv:astro-ph/0411630| . 

[51] U. Alam, V. Sahni, T. D. Saini and A. A. Starobinsky, Mon. Not. Roy. Astron. Soc. 



354, 275 (2004) [astro-ph/0311364] ; U. Alam, V. Sahni and A. A. Starobinsky, JCAP 



0406, 008 (2004) larXiv:astro-ph/0403687l ; JCAP 0702, Oil (2007) | arXiv:astro-ph/0612381 1 ; 



S. Nesseris and L. Perivolaropoulos, JCAP 0701, 018 (2007) | arXi v:astro-ph/0610092] ; 
P. U. Wu and H. W. Yu, Phys. Lett. B 643, 315 (2006) | arXiv:astro-ph/0611507| ; H. K. Jas- 
sal, J. S. Bagla and T. Padmanabhan, Mon. Not. Roy. Astron. Soc. 405, 2639 (2010) 
|arXiv:astro-ph/0601389] . 

[52] Edited by M. Abramowitz and 1. A. Stegun, Handbook of Mathematical Functions with For- 
mulas, Graphs, and Mathematical Tables (1964). 
[53] A. A. Starobinsky, JETP Lett. 68, 757 (1998) [Pisma Zh. Eksp. Teor. Fiz. 68, 721 (1998)] 

28 



|arXiv:astro-ph/981043l] ; D. Huterer and M. S. Turner, Phys. Rev. D 60, 081301 (1999) 
|arXiv:astro-ph/980813"3] ; T. Nakamura and T. Chiba, Mon. Not. Roy. Astron. Soc. 306, 696 
(1999) [arXiv:astro-ph/9810447| . 



[54] V. Sahni and A. Starobinsky, Int. J. Mod. Phys. D 15, 2105 (2006) [arXiv:astro-ph/0610026| . 

[55] A. A. Starobinsky, Grav. Cosmol. 6, 157 (2000) |arXiv:astro-ph/9912054l . 

[56] A. A. Starobinsky, Sov. Astron. Lett. 4, 82 (1978). 

[57] B. Feng, M. Li, Y. S. Piao and X. Zhang, Phys. Lett. B 634, 101 (2006) 
|arXiv:astro-ph/0407432] ; E. Ehzalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 043539 
(2004) [hep-th/0405'034l ; B. Feng, X. L. Wang and X. M. Zhang, Phys. Lett. B 607, 35 (2005) 
|arXiv:astro-ph/0404224] ; X. F. Zhang, H. Li, Y. S. Piao and X. M. Zhang, Mod. Phys. Lett. 
A 21, 231 (2006) [arXiv:astro-ph/0501652| ; M. z. Li, B. Feng and X. m. Zhang, JCAP 0512, 
002 (2005) |arXiv:hep-ph/0503268| ; G. B. Zhao, J. Q. Xia, M. Li, B. Feng and X. Zhang, 
Phys. Rev. D 72, 123515 (2005) |arXiv:astro-ph/0507482] . 

[58] C. Corda, Int. J. Mod. Phys. D 18, 2275 (2009) |arXiv:0905.2"502 ' [gr-ac]]. 



29 



